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The spin-1/2 ferromagnetic-antiferromagnetic alternating Heisenberg chain with ferromagnetic next- 
nearest- neighbour (NNN) interaction is investigated. The ground state is the Haldane phase for weak 
NNN interaction, and is the ferromagnetic phase for weak antiferromagnetic interaction. We find a series 
of topologically distinct spin-gap phases with various magnitudes of edge spins for strong NNN interaction. 
The phase boundaries between these phases are determined on the basis of the DMRG calculation with 
additional spins which compensate the edge spins. It is found that the exact solutions with short range 
antiferromagnetic correlation on the ferromagnetic-nonmagnetic phase boundary is representative of each 
spin gap phase. 

KEYWORDS: frustration, alternating chain, Haldane phase, topological spin-gap phase, edge spin, exact solu- 
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1. Introduction 

Physics of the topological phases of matter is one of the 
most attractive fields of contemporary condensed matter 
physics. Among them, the Haldane phase of the integer- 
spin antiferromagnetic Heisenberg chain^'^ is one of the 
simplest but non-trivial topological phases in quantum 
magnetism. As an extension of the concept of the Hal- 
dane phase, one of the authors introduced the spin-1/2 
Heisenberg chain which has two different exchange in- 
teractions Ja and Jp alternatingly. The ground state of 
this model interpolates those of the spin-1/2 Heisenberg 
chain and spin-1 Heisenberg chain. ^ As long as Ja > Jf 
and Ja is antiferromagnetic, this ground state is adia- 
batically connected to the Haldane phase of the spin-1 
antiferromagnetic Heisenberg chain without closing the 
bulk energy gap. Therefore, this ground state belongs to 
the same topological phase as the Haldane state. 

The quantum magnetism in frustrated spin systems 
is another rapidly developing field of condensed mat- 
ter physics.^'* Although topological phases in two- 
dimensional magnetic systems are expected in frus- 
trated quantum magnets,^' ^'^ the effect of frustration 
on the one-dimensional topological phases has been less 
studied. In the present work, we investigate the ef- 
fect of frustration on the ground state of the spin- 
1/2 ferromagnetic- antiferromagnetic alternating Heisen- 
berg chain. The frustration is introduced by the next- 
nearest-neighbour ferromagnetic interaction Jl . The lat- 
tice structure is shown in Fig. 1. This structure can be 
also regarded as a frustrated ladder^^"-'^'^ with ferromag- 
netic leg Jl, ferromagnetic rung Jp, and antiferromag- 

* E-mail: hida@mail.saitama-u.ac.jp 



netic diagonal interaction Ja. A frustrated ferromagnetic 
ladder compound with this structure has been synthe- 
sized recently.^'' Although the ground state of this mate- 
rial is ferromagnetic, the frustration effect would be en- 
hanced, if the exchange interactions can be modified. As 
a result, the ground state can change into more exotic 
nonmagnetic phases due to the interplay of frustration 
and quantum fiuctuation. Hence, it is worthwhile to in- 
vestigate the possible ground states of this model varying 
the exchange interactions freely. 



,<>-r^ — — — o — o 

Jp^ Va/ \ / \ / \ / \ 

-O-r-O — o — o — o — o 



Fig. 1. Structure of the frustrated ferromagnetic- 
antiferromagnetic alternating Heisenberg chain. 



The main results of the present paper is summarized in 
the phase diagram of Fig. 2. In addition to the ferromag- 
netic (F) phase and the Haldane (H) phase, we find inter- 
mediate phases with various magnitudes of edge spin Se 
in the region with strong ferromagnetic NNN interaction 
Jl. The difference in the magnitude of edge spin reflects 
the topological difference of the bulk ground states. '^"^ A 
well-known example is the edge spins with magnitude 1/2 
in the Haldane phase of the spin-1 Heisenberg chain. 
This state is topologically distinct from the trivial spin 
gap phases such as the dimer phase or the large-D phase 
which has no edge spins. ^■'^ 

This paper is organized as follows. In the next section, 
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Fig. 2. Ground state phase diagram. The ferromagnetic phase, 
the Haldane phase, and the intermediate spin-gap phases with edge 
spin Sc in open chain Hamiltonian (18) arc indicated by F, H, and 
ISo 1 respectively. The stabiUty limit of the ferromagnetic phase (17) 
is shown by the dotted line. The open squares are the boundary 
between diflferent spin-gap phases determined numerically in §4.3. 
The filled circles are the "special points" defined by Dmitriev et 
28, 29 rpj^g solid curves axe guides for the eye. 



the model Hamiltonian is presented. Starting with the 
classical analysis in §3, the ground state phase diagram 
is determined using the numerical diagonalization and 
DMRG technique in §4. The relation to the exact solu- 
tion on the ferromagnetic-nonmagnetic phase boundary 
is also disctissed. The last section is devoted to summary 
and discussion. 

2. Hamiltonian 

We aim to clarify the nature of the ground state of the 
spin- 1/2 Heisenberg chain described by the Hamiltonian 

L 2L 

'Hl = Y^ {JfS 



(1) 



where Si is the spin- 1/2 operator on the ^-th site. In 
the present work, we concentrate ourselves to the case 
Jf, Jl < (ferromagnetic) and Ja > (antiferromag- 
netic). If we assume the periodic boundary condition 
S1+2L = Si, the total number of spins is 2L. In the 
following sections, however, we also consider several dif- 
ferent open chains with additional spins suitable to eluci- 
date the nature of different phases and phase transitions 
between them. 

3. Clcissical Ground-State Phases 

Before analyzing the quantum ground states of Hamil- 
tonian (1), we examine the ground states of its classical 
version for comparison. The classical Hamiltonian "H^' is 
given by replacing the spin 1/2 operator Si by a classical 
vector with length S vo.'Hl- 



S.I Ferromagnetic Phase 

We rewrite T-L'j] in the following complete-square form 

L r 



1=1 



A{{S2l-l-XS2l-{l-X)S2l+lf 
+ {S2l+2 — XS21+I — (1 — X)S2l)^} 

+ Ef, (2) 



+ B{S2i — 821+1)'^ 



-(Jf + 2Jl)/2, B 



where A = Jf/(Jf + 2Jl), A 
-2 Jl Jf/( Jf + 2 Jl) - Ja, and 

Ef = LS^{2Ji^ + Jf + Ja). (3) 

As long as A > and B > 0, the ground state is the 
ferromagnetic state where Si is independent of i, and 
the ground-state energy is Ep. A is always positive for 
Jf > and Jl > 0. Hence the condition B > written 
as 

1 

~2Jl 

gives the ferromagnetic phase. 



^1 1 
~ Jf Ja 



(4) 



3.2 Antiferromagnetic Phase 
We rewrite Hf; in the following complete-square form 



1 ^ 



A{{S2l-l+XS2l-2 + {l-X)S2lf 
+ {S2l + XS2l+l + {l-X)S2l-lf} 



+ B{S2l-2 + S2l-iy 



+ Eaf, 



(5) 



where A = 2Jl/Jf, A = -J|/[2(2Jl - Jf)], B = Ja -|- 
2JlJf/(2Jl - Jf), and 

Eaf = LS^i2A -Jf- Ja)- (6) 

Equations A > and B > mean that the ground 
state is the antiferromagnetic (AF) state where = 
— Si for all I, and the ground state energy is Eaf- These 
conditions reduce to 



1 1 1 

> — H 

2Jl Jf Ja 

for the antiferromagnetic (AF) phase. 



(7) 



3.3 Double-Period Antiferromagnetic Phase 
We rewrite Hf; in the following complete-square form 

L r 



A{{S2l-2+XS2l-l + {l-X)S2lf 

+{S2i-i + XS21-2 + (1 - A)S'2;-3)'} 



+ B{S2l-2 + «S'2i-l)^ 



-Edaf 



(8) 



where A = Jf/(2Jl), A = 2J2/(2Jl - Jf), B = Ja - 
2JlJf/(2Jl — Jf), and the ground-state energy -E^daf is 
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given as 



-E'daf = —LS^{2Ji, — Jf + Ja)- 



(9) 



Equations A > and B > mean that the ground 
state is the antiferromagnetic state with doubled peri- 
odicity (DAF state) where = iS'2;-i = —S21-2 = 
for all I, and the ground state energy is Euaf- 
These conditions reduce to 



1 

'2Jl 



> 



1 _^ 1 
Jf Ja 



(10) 



for the double-period antiferromagnetic (DAF) phase. 

3.4 Spiral Phase 
We rewrite in the following complete-square form 

1 ^ 
1=1 

+{^821 + PS21+1 + aS2i+2?) + Es, (11) 

wit h a = V JaJl/Jf, P = (1/2)VJaJf/Jl, 7 = 
--\/JlJf/Ja, and 

E, = -LS'(:h^ + ^ + ^). (12) 
V Jf 4Jl Ja / 

A ground-state spin configuration is given if all the 
squares in eq. (11) are zero and then the ground state 
energy is £^s- In the ground state, both S21-1 and S21+2 
are on the plane spanned by and S21+1, and then all 
the spins are in the same plane. Setting the plane as the 
a;?/-plane, the spins are given as 



where 



with 



{Scos6n,Ssm6n, 0), 



021 =00 + W, 

021+1 = eo + {i + i)e + <p 



(13) 



(14) 



a 



cos 



6 



i -|- cos 



2 0:7 
-1 



■0" 



2aP 



(15) 



and an arbitrary ^o- This spin configuration shows that 
the system is in the spiral (S) phase with period 9. The 
ground state is a commensurate spiral state ii O/n is a 
rational number, and otherwise an incommensurate one. 
For example, in the case of Jl = — Ja/2 and Jp — — Ja, 
we have a commensurate spiral state with 9 = (j) = -k /2> 
and Es = -(3/2)JaS'2J. 

The spiral spin configuration exists if and only if 
I cos^l < 1 and | cos{9+(j))\ < 1 for eq. (15). These condi- 
tions restrict the values of the exchange parameters and 
determine the phase boundaries between the spiral phase 
and the other phases. Rewriting the conditions by the ex- 
change parameters, we find that the spiral phase is all the 
area except for the F, AF, and DAF phases respectively 




Fig. 3. Classical ground state phase diagram. It includes the fer- 
romagnetic (F), antiferromagnetic (AF), double-period antiferro- 
magnetic (DAF), and spiral (S) phases. 



given by eqs. (4), (7), and (10). The spiral spin configu- 
ration reduces to the F, AF, and DAF configuration on 
the boundaries as expected. 

The classical ground state phase diagram is summa- 
rized in Fig. 3. Comparing this classical phase diagram 
to the quantum one of Fig. 2, the ferromagnetic phase 
remains in the same area. The spiral phase together with 
the AF and DAF phases seem to turn into the Hal- 
dane phase with spin gap and short range order owing to 
quantum fiuctuation. The intermediate phases between 
the ferromagnetic and Haldane phases are proper to the 
quantum system, and will be examined in detail in the 
next section. 

4. Ground State Phase Diagram 

4-1 Strong Coupling Limit : | Jl|, Ja<C| Jf| 

In the strong Jp limit, S21-1 and <S'2i form a triplet pair 
which is described by the spin-1 operator Si = S21-1 + 
821- Then the original Hamiltonian (1) is mapped onto 
an eff'ective spin-1 Heisenberg chain 



(16) 



1=1 



with effective exchange parameter Jeff = (2Jl -I- Ja)/4. 
The ground state is a ferromagnetic state for Jl < 
— Ja/2, while it is a Haldane state for Jl > — Ja/2. 

As long as the spins S'2i-i and form a composite 
spin Si with magnitude 1, the Haldane state has an un- 
ambiguously topological nature. ^"^ Reflecting this fact, 
the edge spins with magnitude 1/2 appear at the open 
edges in the Haldane phase. However, actual spins iS'2;-i 
and ^2; in the Hamiltonian (1) are independent degrees 
of freedom, although they interact with each other. If 
we cut the chain between 52; -1 and ^2;, no edge spins 
appear at the boundary in the Haldane phase. There- 
fore, we can also regard this phase as a trivial spin-gap 
phase. Thus, the definition of the terms "topological" 
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Fig. 4. Jp-dcpcndcncc of the scaled gap LAE with the periodic 
boundary condition for Jl = —1.5. 



and "trivial" is rather arbitrary in the present model. 

To fix the terminology, we choose the pairs of spins on 
the Jf bonds as building blocks of the bulk of our chain 
and call the phase which adiabatically continues to the 
Haldane phase of the spin-1 chain in the limit Jf — >■ — oo 
with Jl < — Ja/2 as "topological". 

4.2 FerrorriMgnetic- Nonmagnetic Phase Boundary 

The stability limit of the ferromagnetic phase can be 
determined from the requirement that the single magnon 
excitation energy vanishes^^ as 

2Ja Jl 



Jf = Jp = ~ 



2 Jl + Ja ■ 



(17) 



This is plotted by a dotted line in Fig. 2. This is the 
same as the classical phase boundary given by (4). It 
is also numerically confirmed that no partial ferromag- 
netic phases^^^^"^ appear between the ferromagnetic and 
nonmagnetic phase s. 

4-3 Intermediate Spin- Gap Phases 

In the following numerical calculation, we fix the en- 
ergy unit by setting Ja = 1 without loss of generality. 
Figure 4 shows the Jp-dependence of the scaled lowest 
single-triplet energy gap LAE for Jl = —1.5 and Ja = 1 
calculated by the nimierical diagonalization method with 
the periodic boundary condition. The size-dependence of 
the scaled gap LAE is weak around Jp ~ —1.3 and LAE 
increases with L on the both sides of this point. This 
suggests that a phase transition between two different 
spin-gap phases takes place around this point. Namely, 
an intermediate spin-gap phase, which is different from 
the Haldane phase, exists between the Haldane and fer- 
romagnetic phases. 

To get more insight into the properties of this inter- 
mediate phase and to determine the phase boundary be- 
tween the intermediate and Haldane phases more pre- 
cisely, we employ the DMRG method for the open chain 



"■"i Jl ^2 Jl 



S2L-2 J S2L 



'^A\-^f/ Ja\Jf' Ja\Jp' J^Jp" J\ 

O-r<>-r-0 O-r-O 

^1 ^3 >32L-1 ' 1 



Fig. 5. Structure of frustrated ferromagnetic-antifcrromagnctic 
alternating chain described by the Hamiltoniaji (19) with spins Ti 
and Ti added on both ends. The total number of spins is 2L + 2. 
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(L+1)AE 



J^=1 Jl=-1.5 



■ L=24 ♦ L=72 

□ L=36 ^ |_=84 

• L=48 

O L=60 



L=96 




Fig. 6. Jp-dependence of the scaled gap with Jl = —1.5 for the 
geometry of Fig. 5. 



Hamiltonian 

L L-l 2L-2 

'HI = ^ JpS2;-iS2; -|- ^ JAS21S21+1 + ^ Ji,SiSi+2 

1=1 1=1 i=l 

(18) 

with the total number 2L of spins. This Hamilto- 
nian naturally continues to the open spin-1 Heiscnberg 
chain in the limit of Jp — >■ —00. Therefore, the quasi- 
degeneracy of the ground state due to the edge spins with 
magnitude 1/2 is expected in the Haldane phase. In the 
present model, however, the lowest magnetic excitation 
gap is extremely small with this boundary condition in 
the whole interval of Jp = — 1.5 < Jp < for Jl = —1.5. 
This suggests that the ground states are quasi-degenerate 
even in the intermediate phase, although the precise anal- 
ysis of the size-dependence of the scaled gap is difficult 
because of the smallness of the energy gap. This is in 
contrast to the conventional Gaussian transition between 
two spin-gap phases such as the Haldanc-dimer transition 
in dimerized spin-1 Heisenberg chains.^'*' 

To resolve the quasi-degeneracy of the ground state 
and the low-lying excited states in the Haldane phase, 
we add two spins Ti and Ti on the both ends of the 
chain as 



HI=HI + Ja{TiSi + S2lTi) 
+ Jl{TiS2 + S2L-1T1). 



(19) 



This geometry is shown in Fig. 5. It should be noted 
that the additional spins Ti and Ti are just added to 
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compensate the edge spins in the Haldane phase to elu- 
cidate the nature of the original open chain W[ which 
tends to the spin-1 open chain in the limit of Jf — >■ — oo. 
In this case, the total number of spins is 2(L + 1). The 
behavior of the scaled gap {L + 1)AE is shown in Fig. 
6. It increases with the system size in the Haldane phase 
as expected, while it still decreases with the system size 
in the intermediate phase. This implies that the edge 
spins in the intermediate phase are not fully compen- 
sated by the added spins Ti and Ti in (19). Actually, 
for the Hamiltonian (19), the lowest energy state with 
i;-component of the total spin S^^^ = 1 and the ground 



state with S'f, 



are quasi-degenerate but the lowest 



energy state with S^g^ = 2 is separated by a gap of 0(1) 
in the intermediate phase. Namely, the ground state is 
quasi-degenerate with the lowest energy state with total 
spin ^tot = 1 as the Kennedy triplet in the open S = 1 
antiferromagnetic Heisenberg chain. We also find that 
the ground state of the chains with no additional spins 
(18) is quasi-degenerate with the lowest energy state with 
5^ot = 2 in the intermediate phase. 

To reveal the spin structure of these states, we plot the 
local expectation values {S^) in the lowest energy state 
of (18) with S'tot = S'tot = 2 in Fig. 7 for Jl = -1.5 
and Jf = —1.4. The accumulated magnetization Mi = 
X]fc=i ("^fc) is ^^^'^ plotted. This plot clearly shows that 
edge spins with magnitude 1 appear on the both ends 
of the chain. At the left (right) end, the magnetization 
is mostly localized on the odd-th (even-th) sites. This 
observation also suggests the possible emergence of the 
phases with larger edge spins for larger |Jl| where the 
odd-th (even-th) spins are strongly correlated with the 
leftmost (rightmost) spin ferromagnetically due to the 
next-nearest-neighbour interaction Jl. 

Motivated by this speculation, we examine the follow- 
ing Hamiltonian 

-Hll, =ni + JaTiSi + JlTiS2 + ^ JlTi-iTi 

1=2 



+ JAfiS2L + JlTi^2L-1 + Yl Ji^Ti.^fi (20) 



Ji_=-1.5 Jf=-1.4 Ja=1 L=72 
■ M21 

□ M2|_1 

1 jl^^^^^^J ^iiii T Trmr 

• <S\|> 

O <S 2|_1> 



1 




50 



Fig. 7. The magnetization profile {Sf') and accumulated mag- 
netization Mi for Jl = —1.5 and Jp = —1.4 . The filled and open 
symbols represent the values for i = 21 and 2/ — 1, respectively. 



(>----0-^>-H> — c> — c>-^, 

(>^K>-p<3 d 0-i-0--:-0 

^1 ^3 ^2L-1 '1 '2 ' L„ 



Fig. 8. Structure of frustrated ferromagnetic-antiferromagnetic 
alternating chain described by the Hamiltonian (20) with spins T; 
and Ti(l = l,...,Le) added on both ends. The total number of 
spins is 2L -f- 2Le . 



(L+2)AE J^=1J^=_1.5 




Fig. 9. Jp-dependence of the scaled gap with Jl = —1.5 for the 
geometry of Fig. 8 with Le = 2. 



with additional 2Le ferromagnetically coupled spins Ti 

and Ti to the both ends of the open chain (18). The lat- 
tice structure is shown in Fig. 8. The total number of 
spins is 2L + 2Le. When a ground state of HI has edge 
spins with magnitude Sc = ic/2, they arc fully compen- 
sated by the added spins in T-LLLe- Hence the phase of 
"H^ is identified by 5e. 

In what follows, we call the intermediate phase which 
has the edge spin with magnitude in (18), the Ig^,- 
phase. The I1/2 phase is the Haldane phase. The Jf- 
dcpendence of the scaled energy gap is shown in Fig. 9 
for Jl = —1.5 and = 2. In this case, the scaled gap 
increases with the system size in the intermediate phase, 
while it decreases with the system size in the Haldane 
phase. This implies that the edge spins in the intermedi- 



ate phase of the Hamiltonian (18) are fully compensated 
by the added spins as expected. This confirms that the 
magnitude Sc of the edge spin is unity in this phase. 

The Jp-dependence of the scaled energy gap is shown 
in Fig. 10 for Jl = —2, and = 1,2 and 3. The same 
plot is shown in Fig. 11 for Jl = —2.5 and Lo = 2,3 
and 4. The scaled energy gap increases with L if the edge 
spins are compensated by the added spins. For Jl = —2.0 
, we find three phases with edge spins = io/2 = 1/2,1 
and 3/2, while for Jl = —2.5 , we find four phases with 
Se = Le/2 = 1/2,1,3/2 and 2 although the data for 
Lc = 1 are not shown. 

Following the scheme of the phenomenological renor- 
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(L+LJAE 




Fig. 10. Jp-dependence of the scaled gap with Jl 
geometry of Fig. 8. Lc = 1, 2 and 3. 



(L+LJAE 




-2. 




-<1-<1— «- 



-<]- 



Jl = -3 -2.5-2-1.75-1.5 -1.3 

L3=1 • ■ ♦ A T ^ 

L.=2 O □ O A V < 



-2 for the 



0.005 



Fig. 11. Jp-dependence of the scaled gap for Jl 
the geometry of Fig. 8.Lc = 2, 3 and 4. 



-2.5 with 



0.01 2/(N,+N2) 



Fig. 12. Extrapolation procedure for the critical values of Jp 
for Jl = —1.5 determined from the intersection point of the scaled 
gap for N = Ni and N = N2 where N = 2L + 2Le is the total 
number of spins. The filled symbols and open symbols are the data 
for Le = 1 and La = 2, respectively. The right-directed triangles 
on the vertical axis are the extrapolated values. 



0.1 

AE 



0.05 



-•Ja-I I, phase oJl= 


-2Jf=-1.1 


I3/2 phase nj^= 


-2.5Jf=-1.15 - 













0.0005 



n2 0.001 



malization group, we determine the finite size phase 
boundary from the intersection point of the scaled gap for 
the Hamihonian (20). To determine the lSa~^s^+i/2 phase 
boundary, we use the energy gap for Hll^ and 'Hll^+i 
with Lc — 25c. After extrapolating to the thermody- 
namic limit L — > 00, the both results coincide within 
2 X 10~^ even in the worst case for the data shown in 
Fig. 2. The extrapolation procedure is shown in Fig. 12 
for the Haldane-Ii phase boundary with Jl = —1.5 . In 
Fig. 2, we plot the lsa~^Sa+i/2 boundary determined from 
^i2So' because its size dependence is weaker than that 
determined from HL2Sa+i empirically. 

The size-dependence of the bulk energy gap in the mid- 
dle of Ii and I3/2 phases are plotted against the system 
size in Fig. 13. It is seen that the magnitude of the energy 
gap is small but finite in these phases. 

4-.4 Valence Bond Structure of the Intermediate Phases 
The structures of the intermediate phases are rather 
enigmatic. For example, the simplest possible structure 
of the Ii-phase speculated from the magnitude of the end 
spins would be the one shown in Fig. 14(a). To examine 
the possibility of this structure we have measured the 
spin correlation functions {S2i+2j+iS2i) with —4 < j < 3 



Fig. 13. Size-dependence of the energy gap for Jl = —2, Jp = 
— 1.1 (Ii-phase) with Lc = 2 and Jl = — 2.5, Jp = —1.15 (I3/2" 
phase) with L^ = 3. 



using DMRG near the center of the chain. 

Contrary to the above intuition, we find that 
{S21S21+1) has the largest amplitude not only in the Hal- 
dane phase but also in the Ii-phase as shown in Fig. 15 
for Jl = —1.5 and Ja = 1. Nevertheless, it is decreas- 
ing with Jp and the correlation {S21S21+3) is increasing 
with Jp. Furthermore, a cusp-like behavior is observed in 
all short range correlation functions at the phase bound- 
ary. Thus, we speculate that the dimer configuration of 
Fig. 14(a) resonate with other longer range valence bond 
structures to lower the energy of the Ii phase. Two ex- 
amples of such dimer configurations are shown in Fig. 
14(b) and (c). They are compatible with the spin gap for 
the Hamiltonian (20) with Lc = 2. 

To confirm this picture, we consider the model with 
additional antiferromagnetic bonds Jad to T-Lll^ with 
Lc = 2 as shown in Fig. 16. The Hamiltonian is given 
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■Jl _^_jJl "Jl "Jl "Jl 





Fig. 14. Structure of the Ii phase with = 2. The ovals corre- 
spond to the valence bond. 



-0.560, 



'0000000000000000000000000^ 

Jl=-1 .5, Ja=1 

O '^S2[_7S2|> O ^S2|S2|+1> 

<S2|_5S2|> • <S2|S2|+3> 

^ <S2|_3S2|> A <S2|S2|+5> 

<S2|_iS2|> ^ <S2|S2|+7> 



-1.5 



-0.5 



Jf 



Fig. 15. Short range spin-spin correlation functions for Jl = 
— 1.5 calculated by the DMRG method. The system size depen- 
dence is negligible except at the H-Ii phase transition point (dotted 
line) where the values extrapolated to L — > 00 are plotted. 



X) — O — ^0 — K) — ^5 — ^ 

Si So T. To 



Fig. 16. Structure of the frustrated ferromagnetic- 
antiferromagnetic alternating chain described by the Hamiltonian 
(21) with The total number of spins is 2L -|- 4. 



by 



L-2 



= 'Hl2 + JaA ^ S21S2I+3 



1=1 



J: 



ad 



T2S1 + T1S3 + T1S2L-2 + T2S 



2L 



(21) 



In the limit of large Jad, it is obviotis that the dimer 
configuration of Fig. 14(a) is the ground state of (21). 
The DMRG result for the Jad-dependence of the energy 



AE 



J^=1 Jl=-1.5 Jp=-1.4 



Fig. 17. Jad-dependence of the energy gap of the Hamiltonian 
(21) with Ja = 1, Jl = —1.5 and Jp = —1.4 calculated by the 
DMRG method. The gap extrapolated to the thermodynamic limit 
L — ^ 00 are shown. 



edge spin Sg 



Ji ®2 Jl ®' 



S2L-4 J S2L_2 J S2L 
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Ja\Ji?* Ja\Jf" , . — J 

O o 6 h h b edge spin Se+1 

S, S2 S3 Jl Sj Jl S2L-3 



Fig. 18. Local modification of "H? . 



gap is shown in Fig. 17. This shows that the ground state 
for large Jad adiabatically continue to that with Jad = 
without closing the energy gap. 

The distinction between spin-gap phases with differ- 
ent edge spin is ascribed to the topological difference in 
the bulk ground states of (18)."^"^'^^ According to rcf. 
27, only two topologically distinct spin-gap phases are 
possible in the spin- 1/2 chains without translational in- 
variance. 

To see this explicitly for the present model, we modify 
the Hamiltonian "HI into the Hamiltonian T-C^°^ given 

by, 



L-l 



L-2 



mod 



JfS21-\S21 + JAS21S21+I 



1=2 

2L-4 



1=2 



J'LSiSi+2 + JhSiS2 + Jl'S'2<S'3 



i=3 



+ JLS2L-2S2L-I + JlS2L-iS2L- 



(22) 



as shown in Fig. 18. This modification consists of the 
following local change of the exchange constants: The 
exchange constants of S1-S2 and S2L-1-S2L bonds are 
changed from Jp to Jl. Those of S2-S3 and S2L-2- 
S2L-1 bonds are changed from Ja to Jl. Those of Si- 
S3, S2-S4, S2L-2-S2L, and S2L-3-S2L-1 bonds are set 
equal to 0. If the ground state of HI has edge spin s with 
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Fig. 19. The accumulated magnetization Mi on the magnetic- 
nonmagnetic phase boundary (17) at tiie special points Jl/Ja = 
— m/2('m = 3,4 and 5) in the ground state witii Stot = S^^^ = 
m — 1. The filled and open symbols represent the values for i = 21 
and 21 — 1, respectively. 



magnitude Sc, that of the modified Hamiltonian H™°'^ 
has edge spin s with magnitude 5*6 + 1 for ferromagnetic 
Jl. This implies that /5,-phase and /g^+i-phase are ac- 
tually the same phase , since this modification only con- 
cerns the local bonds and the bulk ground state should 
remain unchanged, as long as it is unique and gapful ex- 
cept for the quasi-dcgeneracy due to the edge spins. In 
this context, all the ground states with integer Se form 
a single phase and those with half-odd-integer Se form 
another single phase. However, as discussed in §4.1, it is 
arbitrary which of the two phases should be called "topo- 
logical" and which one "trivial" . Therefore, we regard the 
Is^ phase with half-odd-integer 5*0 which is connected to 
the Haldane phase of the spin-1 chain "topological" and 
that with integer "trivial" phase. Within the parame- 
ter space of the present model, however, we cannot move 
from the Ig^ -phase to the Is^+i-phase without passing 
through the l5^_|_i/2-phase. Therefore, all Is^-phases with 
different Sc are separate phases in the present model. 
4-5 Edge Spin State on the Ferromagnetic- 
Nonmagnetic Phase Boundary 
The exact ground state on the ferromagnetic- 
nonmagnetic phase boundary (17) is obtained by 
Dmitriev et al.^^'^^ for the Hamiltonian 

Jf 

"H/f = -—{S1S2 + S2L-1S2L) + JFS21-1S21 
^ 1=2 

L-l 2L-2 
+ ^ JAS21S21+I + ^ Ji,SiSi+2, (23) 
1=1 j=l 

where the exchange couplings of the endmost bonds are 

halved to allow the exact solution. The ground state 
for J^/Ja = —to/2, where m is a positive integer, has 
a short range antiferromagnetic order. Therefore, it is 
expected that this solution continues to the spin-gap 
ground state. These points are called special points. On 
the contrary, the ground state for Jl/Ja 7^ — m/2 has a 
long range spiral spin correlation with wave length equal 
to the system size, which would correspond to the ferro- 



magnetic state with a single twist. From the continuity 
consideration, it is plausible that both types of eigen- 
states are actually present for arbitrary values of Jl/Ja 
on the ferromagnetic-nonmagnetic phase boundary and 
they are degenerate in the thermodynamic limit. The 
presence of the macroscopic degeneracy is proven on the 
special points, while it is suggested by the numerical re- 
sults for other points on the ferromagnetic-nonmagnetic 
phase boundary.^^ Thus, we examine the solutions with 
short range antiferromagnetic order at the special points 
Jl/Ja = -m/2. 

Figure 19 shows the accumulated magnetization Mj 
for total spin Stot = St^t = m — 1 calculated using the 
method of ref. 30. This shows that these solutions have 
edge spins with S^ = (m — l)/2 on the both ends of the 
chain. As shown in the phase diagram of Fig. 2, each 
special point with m = 2Se + 1 belongs to the /s^, -phase. 
Thus, we may regard that the special point solution with 
m is representative of Is^ phase for each 5c = (m — l)/2. 

Among them, the ground state solution on the special 
point with m = 2 is an isolated dimer state on the Ja 
bonds. This state remains the ground state all the way 
along the line Jp = 2 Jl for Jl > — Ja- This line belongs 
to the Haldane phase in the phase diagram Fig. 2. This 
confirms that the special point solution with to = 2 is 
representative of the valence bond structure of the Hal- 
dane phase. 

5. Summary and Discussion 

Ground-state phases of the spin-1/2 ferromagnetic- 
antifcrromagnctic alternating Hciscnbcrg chain with fer- 
romagnetic NNN interaction arc investigated. 

In addition to the conventional Haldane phase and the 
ferromagnetic phase, we confirmed the presence of in- 
termediate spin-gap phases which has edge spins with 
various magnitudes on both ends. These phases are clas- 
sified into two topologically distinct phases according to 
whether the edge spin 5c is an integer or half-odd-integer 
in the open-chain Hamiltonian (18). The physical pic- 
ture of each phase is discussed based on the appropri- 
ate modification of the Hamiltonian. The relation to the 
exact solution on the ferromagnetic-nonmagnetic phase 
boundary is also discussed. It is found that each special 
point solution^^' with Jl/Ja = —m/2 has edge spin 
(m — l)/2 and is representative of the Jg^ -phase with 
5e = (to - l)/2. 

Although we have characterized the nature of the spin- 
gap phases in our model as a series of topological phases 
by detecting the edge spins, several types of exotic phases 
are known to exist near the ferromagnetic-nonmagnetic 
phase boundary in other frustrated quantum spin chains. 
Hence, it is worthwhile to examine the possibility of 
their realization in the intermediate phases of the present 
model. 

One of the candidates is the spin nematic phase which 
is a condensed state of bound multimagnons.^^'^^'^^'^^ 
In this case, the lowest bulk excitation should have total 
spin 2 or larger. We have checked numerically that this 
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is not the case in the intermediate phases of the present 
model. 

Spontaneously polymerized phases with broken trans- 
lational symmetry are also candidates.^^'^'^'"^^''^'^"'^''' In 
the open chain, however, one of the symmetry-broken 
state should be automatically selected, if the system size 
is an integer multiple of the periodicity of the ground 
state. Hence, a spatial modulation of physical quanti- 
ties should be observed in the ground state calculated by 
DMRG. Such a spatial modulation is not observed in the 
intermediate phases of the present model. 

The spiral (quasi-)long range ordered phase is also ex- 
cluded because of the presence of the bulk spin gap. Nev- 
ertheless, the short range spiral order might be possible 
even in the spin-gap phases. In this case, the spin gap 
should open at a finite wave number. Although we 
found no evidence for this behavior within the system 
size accessible by the numerical exact diagonalization 
method, the possibility of the spiral short range order 
with the pitch longer than the system size remains pos- 
sible. Even in this case, the ground state phase transition 
does not take place as long as the spiral order remains 
short ranged. 

The physical origin of the edge spin with general values 
of Se may be understood in the following way: For large 
ferromagnetic Jl, the spins on the odd-th (even-th) sites 
are strongly correlated ferromagnetically each other. In 
the nonmagnetic phase, however, this correlation does 
not extend over the whole chain, but is cut into finite 
ferromagnetic clusters fluctuating in position and length. 
The clusters on the even-th sites and odd-th sites are 
correlated antiferromagnetically forming a nonmagnetic 
ground state as a whole. On both ends of the chain, the 
endmost clusters arc pinned to the edges resulting in the 
edge states with spin 5*0. 

The numerical diagonalization program is based on the 
package TITPACK ver.2 coded by H. Nishimori. The nu- 
merical computation in this work has been carried out 
using the facilities of the Supercomputer Center, Insti- 
tute for Solid State Physics, University of Tokyo and Su- 
percomputing Division, Information Technology Center, 
University of Tokyo, and Yukawa Institute Computer Fa- 
cility in Kyoto University. This work is supported by a 
Grant-in-Aid for Scientific Research (C) (21540379) from 
Japan Society for the Promotion of Science. 
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